DISCRETE FLAT SURFACES AND LINEAR WEINGARTEN SURFACES 

IN HYPERBOLIC 3-SPACE 

TIM HOFFMANN, WAYNE ROSSMAN, TAKESHI SASAKI, AND MASAAKI YOSHIDA 

Abstract. We define discrete flat surfaces in hyperbolic 3-space from the perspective 
of discrete integrable systems and prove properties that justify the definition. We show how 
these surfaces correspond to previously defined discrete constant mean curvature 1 surfaces 
m H^, and we also describe discrete focal surfaces (discrete caustics) that can be used to 
define singularities on discrete flat surfaces. Along the way, we also examine discrete linear 
Weingarten surfaces of Bryant type in , and consider an example of a discrete flat surface 
related to the Airy equation that exhibits swallowtail singularities and a Stokes phenomenon. 
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2 tim hoffmann, wayne rossman, takeshi sasaki, and masaaki yoshida 

1. Introduction 

There is the classical Weierstrass representation for minimal surfaces in Euclidean 3-space 
M'', which gives a local conformal parametrization for any minimal surface. It involves choosing 
two holomorphic functions (or perhaps meromorphic functions when considering the surfaces 
more globally) on a Riemann surface. If one restricts to isothermic parametrizations, that is, 
conformal parametrizations that are also curvature lines coordinate systems, then the repre- 
sentation reduces to the choice of just one holomorphic function. Since every minimal surface 
has local isothermic coordinates (away from umbilics), this reduction does not involve any loss 
of generality beyond avoiding umbilic points. 

Once one restricts to isothermic parametrizations, it becomes possible to give a definition 
for discrete analogs of minimal surfaces [6]. These analogs are "discrete isothermic" meshes 
(a definition of this is given later in this paper). They are comprized of planar quadrilaterals, 
which in particular have concircular vertices. 

By a transformation called the Lawson correspondence or T-transformation or Calapso trans- 
formation [14j . one can produce all constant mean curvature (CMC) 1 surfaces in hyperbolic 
3-space H'^ from minimal surfaces in R^. There is a corresponding holomorphic representation 
for those surfaces as well, first given by Bryant [TU] ■ Correspondingly, without loss of generality 
beyond avoiding umbilics, one can restrict to isothermic coordinates in this case also, and one 
has a discrete analog of CMC 1 surfaces in H^, first found by Hertrich-Jeromin [13]. 

In the case of smooth surfaces there is also a holomorphic representation for flat (i.e. in- 
trinsic curvature zero) surfaces in H'^ [11] and this also ties in to the above-mentioned Bryant 
representation, as there are deformations from CMC 1 surfaces in to flat surfaces via a 
family of linear Weingarten surfaces in [T^]. These do not include all linear Weingarten 
surfaces, but rather a certain special subclass called linear Weingarten surfaces of Bryant type 
[l2] [20] . so named because they have Bryant-type representations. 

Thus it is natural to wonder if flat surfaces also have a discrete analog, and we will see here 
that they do. Once this discrete analog is found, a new question about "singularities on discrete 
flat surfaces" naturally presents itself, in this way: Unlike the smooth isothermic minimal 
surfaces in and CMC 1 surfaces in H'^, smooth flat fronts have certain types of singularities, 
such as cuspidal edges and swallowtails (in fact, indirectly, this is what the naming "fronts" - 
rather than "surfaces" - indicates). The means for recognizing where the singularities are on 
smooth flat fronts are clear, and one can make classifications of those surfaces' most generic 
types of singularities just from looking at the choices of holomorphic functions used in their 
representation [17] . However, in the case of discrete flat surfaces, it is not apriori clear where 
the singularities are, nor even what such a notion would mean. Since one does not have first 
and second fundamental forms at one's disposal in the discrete case, one must find an alternate 
way of defining singularities. We aim towards this by defining and using a discrete analog of 
caustics, also called focal surfaces, for smooth flat fronts. For a smooth flat front, the caustic is 
the surface comprised of all the singular points on all parallel surfaces of that flat front. (The 
parallel surfaces are also flat.) Thus the singular set of the flat front can be retrieved by taking 
its intersection with its caustic. In the case of a smooth flat front, the caustic is again a flat 
surface, but this will not quite be the case for discrete flat surfaces. 

We will also present a number of examples of these discrete flat surfaces. In addition to the 
rather simple examples of discrete cylinders and discrete surfaces of revolution, we will also 
discuss a discrete flat surface based on the Airy equation. This example exhibits swallowtail 
singularities and a Stokes phenomenon, similar to that of the analogous surface in the smooth 
case, as shown by two of the authors in [27]. This last example hints at existence of a robust 
collection of discrete flat surfaces with interesting geometric properties yet to be explored. 

Thus, the purpose of this paper is to: 

(1) provide a deflnition for discrete flat surfaces and discrete linear Weingarten surfaces of 
Bryant type in hyperbolic 3-space H'^; 
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(2) give properties of these surfaces that justify our choice of definitions (in particular, as 
smooth flat fronts have extrinsic curvature 1, we identify notions of discrete extrinsic 
curvature of discrete flat surfaces which do indeed attain the value 1); 

(3) show that these surfaces have concircular quadrilaterals; 

(4) study examples of these surfaces, and in particular look at swallowtail singularities and 
global properties of an example related to the Airy equation; 

(5) give a definition of discrete caustics for discrete flat surfaces; 

(6) show that the caustics also have concircular quadrilaterals and that they provide a 
means for identifying a notion of singularities on discrete flat surfaces. 

In Section [2] we describe smooth and discrete minimal surfaces in Euclidean 3-space M'^, to help 
motivate later definitions, and we also give the definition of a discrete holomorphic function, 
which will be essential to everything that follows. In Section [3] we describe smooth CMC 
1 surfaces, and flat surfaces and linear Weingarten surfaces of Bryant type in H'^, again as 
motivational material for the definitions of the corresponding discrete surfaces in Section 31 
We prove in Section |4] that discrete flat surfaces and linear Weingarten surfaces of Bryant type 
have concircular quadrilaterals. Also, Section [4] provides a natural representation for discrete 
flat surfaces which gives the mapping of the surfaces as products of 2 by 2 matrices times 
their conjugate transposes, and we show that this representation applies to the case of discrete 
CMC 1 surfaces as well. The definition for discrete CMC 1 surfaces is already known [T3], but 
the representation here for those surfaces is new. In Section [5] we look at a specific discrete 
example whose smooth analog is equivalent to solutions of the Airy equation, and we look at 
the asymptotic behavior of that surface, which exhibits swallowtail singularities and a Stokes 
phenomenon. In Section [6l we look at normal lines to discrete flat surfaces. With this we 
can do several things. For example, we look at parallel surfaces (which are also discrete flat) 
and show that the area of corresponding quadrilaterals of the normal map equals the area of 
the quadrilaterals of the surface itself, as should be expected, since in some sense the extrinsic 
curvature of the surface is identically equal to 1 (note that the analogous statement is true for 
smooth surfaces with extrinsic curvature 1, infinitesimally) . Then, using distances from the 
surface's vertices to the intersection points of the normal lines, we consider a discrete analog 
of the extrinsic curvature and see that it is 1 in the discrete case as well. Furthermore, those 
intersections give us a means to define discrete caustics, and, as mentioned above, we use those 
caustics to study the nature of "singularities" on discrete flat surfaces, in the final Section [T] 

Acknovirledgements. The authors thank Udo Hertrich-Jeromin for fruitful discussions and 
valuable comments. 

2. Smooth and discrete minimal surfaces in 

A useful choice of coordinates for a surface is isothermic coordinates. Not all surfaces have 
such coordinates, but CMC surfaces in space forms such as M.^ and do have them, away 
from umbilic points. Isothermic coordinates will be of central importance in this paper. 

Another useful tool in the study of surfaces in space forms is the Hopf differential, which 
is defined as Q = {fzz, N)dz'^ , where the surface / is a map from points z in a portion of the 
complex plane C, (•, •) is the bilinear extension of the metric for the ambient space form to 
complex vectors, and N is the unit normal to the surface. When the coordinate z is conformal 
and the surface is CMC, then Q will be holomorphic in z. Umbilic points of the surface occur 
precisely at the zeros of the Hopf differential. 

2.1. The Weierstrass representation for smooth minimal surfaces. Locally, away from 
umbilics, we can always take a smooth minimal immersion / = /(x, y) into M'^ to have isothermic 
coordinates {x, y) in a domain of R^. Let N denote the unit normal vector to /. Then, setting 
z — X -\- iy, the Hopf differential becomes Q — rdz^ for some real constant r, and rescaling the 
coordinate z, we may assume r — I. 



4 



TIM HOFFMANN, WAYNE ROSSMAN, TAKESHI SASAKI, AND MASAAKI YOSHIDA 



Let g be the stereographic projection of the Gauss map N to the complex plane, and set 
g' = dg/dz. As we are only concerned with the local behavior of the surface, and we are allowed 
to replace the surface with any rigid motion of it, we may ignore the possibility that g has poles 
or other singularities, and so the map g : C — > C is holomorphic. Because we avoid umbilic 
points of /, we also know that g' is never zero. Thus the Weierstrass representation is (with 
regarded as lying in the complex plane C) 

f^Re r\2g,l-g',V^{l+g^))cj, cj ^ 9. ^ ^ , 
Jzo dg g' 

Associating (1,0,0), (0, 1,0) and (0,0, 1) with the quaternions i, j and k, respectively, we have 

(2-1) fx^{i~ 9j)j — {i - gj) , fy^{i~ 93)j — {i - 9i) ■ 

9x 9y 

We have converted to a formulation using quaternions here, because this type of formulation 
has been used to define discrete minimal surfaces in M'^ and discrete CMC 1 surfaces in H'^, 
and we wish to make comparisons to those formulations. 

Note that by restricting to isothermic coordinates, we can then determine minimal surfaces 
by choosing just one holomorphic function g. 

2.2. Discrete holomorphic functions. To define discrete minimal surfaces, we use discrete 
holomorphic functions g = gm,n : D <C, where D is the square integer lattice Z^, or a 
subdomain of it. Discrete holomorphic functions are defined as follows: defining the cross ratio 
of 9 to be 

Crjn.n ^ (ffrjj.+ l, ri 9'm,n)i9m+l,n+l 9rn+l,n) {9rn,n+l 9m+l,n+l 

we say that g is discrete holomorphic if there exists a discrete mapping a to M such that 
(2.2) cr„,„ = < , 

'^(m,n) (m,n+l) 

with a(m,„)(m+i,n) = a(m,n+i)(m+i,n+i) and a(^m,n){m,n+i) = a(m+i,n)(m+i,n+i) ioT all quadri- 
laterals (squares with edge length 1 and vertices in D.) See [7]- We call the discrete map a a 
cross ratio factorizing function for g. 

Note that a is defined on edges of D, not vertices. Note also that a is symmetric, that is, 

^(m,n,)(m+l,n) *-*^(m+l,n)(m,n) and Q^(m,n,)(m,n+1) '^(m,n+l)(m,n) ■ 

There is a freedom of a single real factor in the choice of these a(^jyi,n){m+i,n) and ck(m,n)(m,n+i) j 
since we could replace all of them with Xce(m,n){m+i,n) and Xa(^m,n)(m,n+i) for any nonzero real 
constant A, and all relevant properties would still hold. Throughout this paper we use A to 
denote that free factor. 

In the above definition of the cross ratio, we have a product of four terms. Since gm.n G 
these terms all commute, and so we could have written this cross ratio simply as a product of 
two fractions. However, when we later consider the cross ratio for quaternionic-valued objects 
or matrix-valued objects, commutativity no longer holds and the order of the product in the 
cross ratio becomes vital. So, for later reference, we have chosen to write the cross ratio in the 
somewhat cumbersome way above. 

The definition above for discrete holomorphic functions is in the "broad" sense. The defini- 
tion in the "narrow" sense would be that cr^,™ is identically —1 on D (see [5], [?])■ Furthermore, 
note that, unlike the case of smooth holomorphic functions, the discrete derivative or discrete 
integral of a discrete holomorphic function is generally not another discrete holomorphic func- 
tion. 

Let us exhibit some examples of discrete holomorphic functions: 

(1) Let = 1? = {(to, n) | m, n e Z}, and set gm.n = c{m + in) for c a complex constant. 

(2) Let _D = Z^, and set gm.n — e^^™"*"™-* for c a real or pure imaginary constant. One 
could also take the function e'^i'^+'^^n fQj. choices of real constants ci and C2 so that 
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the cross ratio is identically —1, giving a discrete holomorphic function in the narrow 
sense. 

(3) In Section [5] we will describe a discrete flat surface based on a discrete version of 
the power function g = z'^ (7 G R), which we define here. This function is discrete 
holomorphic in the narrow sense. It is defined by the recursion 

(2.3) 7 • g,„,„ = 2m h 2n . 

.9m+l,n ^/m— l,n ^/m,n+l i7m,n— 1 

We start with D — {(to, n) | m, n > 0}. For 7 e (0, 2), the initial conditions should be 

50,0 = , 91,0 = 1 , 50,1 = • 

We can then use (|2.3|) to propagate along the positive axes {9m,o} and {50,™} with 
TO > 1 and n > 1, respectively. We can then compute general gm,n (for both to > and 
n > 0) by using that the cross ratio is always —1. The 9m.n will then automatically 
satisfy the recursion relation (|2.3p . This definition of the discrete power function can be 
found in Bobenko [2]. (It is also found in a recently published textbook [9].) Agafonov 
[1] showed that these discrete power functions are embedded in wedges (see Figure [T]) , 
and are Schramm circle packings (see [30]). Note that, for to G Z and m > 1, 

(rt A\ V2/m -7 

(^■4) 32m, — -7 ^ , ff2m+l,0 — "7 ^^T , 50,n — * ffn,0 , 

\2)m+l V 2jm+l 

where (a)™ = a{a + l)...(a + to — 1) denotes the Pochhammer symbol, and a closed 
expression for general 9m,n is still unknown. We explore this difference equation ()2.3p 
in more detail in Appendix [5] at the end of this paper. 




Figure 1. Domain (left) and image (right) for the discrete power function 



,4/3 



2.3. Discrete minimal surfaces. The representation (|2.ip for smooth minimal surfaces above 
suggests that the definition for discrete minimal surfaces is (see |13) ) 

(2-5) fq- fp = {i- 9pj)j [i ~ 9qj) , 

9q 9p 

where g : D Cis a discrete holomorphic function with cross ratio factorizing function a, and 
p and q are either (m,n) and (n + l,n), or (to, n) and (rn^n +1). This defines the surface / 
up to translations of M?" . The freedom of scaling of a leads to homotheties of /. 

As in the smooth case, where we avoided umbilics, and thus 9' was never zero, we will make 
the following assumption throughout this paper: 

Assumption: gq ~ gp ^ ■ 

Example 2.1. The discrete holomorphic function c(to + in) for c a complex constant will 
produce a minimal surface called a discrete Enneper surface, and graphics for this surface can 
be seen in [6]. 



6 



TIM HOFFMANN, WAYNE ROSSMAN, TAKESHI SASAKI, AND MASAAKI YOSHIDA 



Example 2.2. The discrete holomorphic function e^i™^*^^" for choices of constants ci and 
C2 so that the cross ratio is identicaUy —1 will produce a minimal surface called a discrete 
catenoid, and graphics for this surface also can be seen in [6]. 

3. Smooth CMC 1 surfaces, flat fronts, and linear Weingarten surfaces in 

3.1. Smooth CMC 1 surfaces. Similarly to the case of minimal surfaces, we can describe 
smooth and discrete CMC 1 surfaces in H^. Hyperbolic 3-space H'^, considered in Minkowski 
4-space K.'^^"'^ = {{xq, xi,X2,X3) \ Xj S M} (with Minkowski metric —dxQ + dx\ + + dx^), is 

= {(a;o,a::i,a;2,a;3) e R^'^ | xq > 0, - - X2 - = 1} « 



X = 



tv{X) > 0,dctX = 1 ^ = {F-F^|Fe SL2{C)} , 



Xq + ^3 2:1+ 1X2 
^Xi — iX2 Xq — X3 

where the superscript T denotes transposition. 

A smooth isothermically-parametrized CMC 1 surface (away from umbilic points), has the 
Bryant equation [10 



(3.1) dF^Fl^^ ^9} J' P^SL2{C), 

where 5 is a holomorphic function with nonzero derivative, and the surface is then 

fi = F ■ F'^ . 

3.2. Smooth flat fronts. Starting with a smooth CMC 1 surface /i with lift F as above, 
define 



2\ dz 



1 g 

1 



(3.2) E = F- 
A flat front is then given by 

with unit normal vector field 

(3.3) N = E-(\ \\-E^ 



We know this surface /o is flat, because (see [12]) 

(3.4) dE = E(^J^_, ^'^dz. 

This surface with singularities is actually a front, because jg'p + Iff'l^^ > 0, which means that 
the associated Sasakian metric is positive definite. (See Theorem 2.9 in [21 .) The notion of 
fronts is important in the study of smooth flat surfaces with singularities in H'^. For example, 
it is a necessary notion for considering the caustic of a flat surface with singularities. However, 
it will not play such a direct role in our considerations on discrete surfaces here. Thus, from 
here on out, we will simply consider smooth flat surfaces with singularities, and sometimes will 
even just call them flat surfaces even though they might have singularities. When the front 
property is actually playing a role, we shall parenthetically refer to the word "front". For more 
information about flat fronts, see [17], [18] and [2T] . 

Remark 3.1. Because the off-diagonal terms g' and {g')~^ in (|3.4p are inverse to each other, 
the conditions in [TT for having singular points, cuspidal edges and swallowtails on /o simplify 
to this: 

(1) singular points occur precisely at points where \g'\ = 1; 

(2) a singular point is a cuspidal edge if and only if Im(^) 7^ ("Re" and "Im" denote 
the real and imaginary parts, respectively); 

(3) a singular point is a swallowtail if and only if g" ^ and Im(^) = and R.e((^)') ^ 0. 
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The condition that Im(^) — holds at some point along a singular curve |g'| = 1 is equivalent 
to the curve having a vertical tangent line at that point. 

3.3. The hyperbolic Schwarz map and a special coordinate w. For later reference, we 
can take a new coordinate w such that 

dw — — dz , 

g' 

which is locally well-defined, and still conformal, although not necessarily isothermic. This 
gives 

(3.5) E-^dE = (J 2 ) dw , 



with 



w{z) 

The reason for changing variables from z to w is that Equation (|3.5p now becomes 
(3.6) £,u-qiw)-u = 

with /o — EE^ the hyperbolic Schwarz map (see [27]), where 



(3.7) E : 



"1 



with functions ui, U2 that are linearly independent solutions of Equation (|3.6p chosen so that the 
constant Ui-^{u2) — U2^(ui) will be 1. This equation p.6p with q{w) — w is the well-known 
Airy equation. 



Remark 3.2. Using q and w, the conditions in Remark 13.11 become: 

(1) Singular points: \q\ — 1, 

(2) Cuspidal edge points: Im(^72) 7^ 0, 

(3) Swallowtail points: g„ ^ 0, Im(^) = and Re(^2=^iM.) ^ 0. 

3.4. Smooth linear Weingarten surfaces of Bryant type. We will now give a deformation 
through linear Weingarten surfaces between the surfaces /i and /o described in Sections 13.11 
and 13.21 This deformation was first introduced in [12j . There are numerous ways to choose 
the deformation, and no one way is geometrically more canonical than any other. We will soon 
come back to this issue (Section 13. 5p . However, for now we will simply fix one choice for the 
deformation - the one that deforms the /i and /o as given in Sections 13.11 and 13.21 into each 
other, in accordance with the notations of previous papers ([H], [H], [12], [H])- 

This particular choice will suffice when we switch to investigating discrete surfaces later. 
In fact, the resulting classes of discrete flat surfaces and discrete linear Weingarten surfaces 
of Bryant type, although defined in terms of the deformation, do not actually depend on the 
choice of deformation. (We say more about this in Section |4^ ) Hence the theorems we prove 
about these surfaces also are independent of the choice of deformation. 

We shall refer to this choice of deformation as the "first Weingarten family" (of either /i or 
/o), and the procedure that we follow for constructing it is as follows. Following [T^ . 

• we convert the E in [12], actually notated as "g" there, to (i?-^)^^, 

• then changing the holomorphic function h in [12 to —g for the function g given here, 

• and allowing —oj in [12j to become dz/g' for the function g here, 

• and also changing / and N to (/"^)~^ and (iV^)~^, respectively, 
a linear Weingarten surface of Bryant type in satisfying (see also j20) ) 

(3.8) 2t{H ~l) + {l^t)K , 
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where H and K are the mean and intrinsic curvatures, respectively, is 

ft - {EL){ELf . 

Here E satisfies p.4p and 



All linear Weingarten surfaces satisfying (|3.8p (i.e. of Bryant type) can be consructed in this 
way (see [l2]). 

When t = 0, we use the frame E = E ■ (L|t=o) = F ■ ^ ^ , giving a flat surface /o, as in 
Section [321 When t = 1, we use the frame F ^ E ■ {L\t=i) ^ E ■ '^^Y giving a CMC 1 



surface /i, as in Section 13.11 Thus we have a deformation through linear Weingarten surfaces 
in H^, from CMC 1 surfaces in to flat surfaces in H^. (See also [20].) 

3.5. Geometric non-uniqueness of the deformation in [12]. We now explain in more 
detail why there is non-uniqueness for the choice of deformation through linear Weingarten 
surfaces of Bryant type. The deformation between smooth flat surfaces and smooth CMC 1 
surfaces through linear Weingarten surfaces, given in Section 13. 4i is not uniquely determined 
in any geometric sense, because of ambiguities in the choice of Weierstrass data. We illustrate 
this with two lemmas, both of which are easily verified: 

Lemma 3.3. Given a smooth isothermically-parametrized CMC 1 surface fi in with lift F 
and Weierstrass data g, the transformation 

F^F-B , B = ( P_ Vl eSU2 

will not change the resulting surface fi — F ■ F^ — FB ■ FB , and will change the Weierstrass 
data by 

- pg - q 

9^9 = — ■ 

q9 + p 

Remark 3.4. By Lemma [3.31 different choices for B do not affect /i. However, when B is not 
diagonal, the transformation in the above lemma generally will result in a different deformation 
ft through linear Weingarten surfaces for t < 1, and also in a different flat surface /q. 



Lemma 3.5. Given a flat surface fo with Weierstrass data g' and lift E as in Equation (j3.4p . 
then the transformation 

g-> g + a , 

followed by the transformation 

E^ E-B , 

where B G SU2 is either diagonal or off-diagonal and a is any complex constant, will not change 
the resulting surface fo^E- E^ — EB ■ EB . 

Remark 3.6. Under the transformations of i? and g in Lemma 15751 /o is not affected. However, 
when B is off-diagonal or a is not zero, these transformations generally will result in a different 
deformation ft through linear Weingarten surfaces for t > 0. In particular, the CMC 1 surface 
/i generally will change. 

The two lemmas and two remarks above show that the " first Weingarten family" deformation 
will change when different Weierstrass data is used, even when the original surface under 
consideration does not change. 

To demonstrate that other choices actually do give different deformations, we will also con- 
sider a deformation we call the "second Weingarten family" , given by starting with a flat surface 
fo with given lift E, and then using the lift EB for an off-diagonal B S SU2 instead to make 
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the deformation through hnear Weingarten surfaces of Bryant type. (Any choice of off-diagonal 
B G SU2 will result in the same deformation.) Regardless of whether E or EB is used, we 
have the same surface /o, but E and EB give opposite orientations for the normal vector to 
/o, and we are interested in this particular choice for a second deformation precisely because 
of this orientation-reversing property. The frames E and EB give different families of linear 
Weingarten surfaces when t > 0. Such different deformations of surfaces can be seen in Fig- 
ures [5] (first Weingarten family) and [3] (second Weingarten family), and also in Figures [5] (first 
Weingarten family) and [6] (second Weingarten family) . 

3.6. The deformations with respect to the coordinate w. The first Weingarten family 
is 

/i^) = {EL)(ELf , 

where L is as in Section 13.41 and E solves Equation p.4p . In terms of the new coordinate w 
given in Section [231 since g — J qdw, L takes the form 



(2) 

A second Weingarten family can be given by taking 



B 



1 

-1 




and then 
where 

,0 y V l + t^hh ' J g' 

and E again solves p.4p . Note that there is freedom of choice of additive constant in the 
definition of h, and different constants will give different deformations. 

In terms of the new coordinate w satisfying dw = ydz (and q = {g')^) given in Section [3.31 

(2) 

a second Weingarten family ' has a particularly nice expression for its singular set: 

(2) 

Lemma 3.7. The second Weingarten family fl taken by choosing h = ~w is singular along 
the curve 

\q\^{l + t\w\')^-{l-t)' = 0. 

(2) 

Note that, in particular, could be singular only at points where q — 0. 



Proof. Note that 



{EB)-'d{EB) = -{^^ l]dw 



and 



{2)_fxo + X3 Xi+iX2\ _ I^^^H- -tw \ sTi^T 



\ X^ — iXo Xn ~ X-i / \ 1 ^+lo,,|2 ' 



yXi-tX2 xo~X3j \ ^tw l + t\wY 

where the Xj are now considered as real- valued functions of w and w (and also of t) . Substituting 
'^^fe ~ ^w'^™ ^ iM"*^^ Minkowski norm dx^ + dx2 + dx^ — o?Xq, we find that 

^dw^ + %dw^ + 2^dwdw = dxf + dxj + dx^ - dxl 

has discriminant 

»2 ~ 2t2t = Wi,{u2) - U2^{u,t (|gp(l + t\w\^ f - (1 - tff (1 + t\w\'')-' > , 

where ui and U2 are as in Section [3. 31 Since mi^(u2) ~ U2-^{ui) = 1, the proof is completed. 

□ 
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3.7. Examples. We now give some examples. 

Example 3.8. Take any constant q G C \ {0} in Equation (|3.5p . Then we can take g = ^fqz = 
qw, and so both the coordinates z and w will be isothermic if (7 S M, and we now assume q is 
a positive real. Let us use the coordinate w, and then E as in Section [3.31 can be taken as 



F = q- 



^ ^ 1/4 / cosh(^w) ^sinh(^u;)\ 
\^sinh(y/gw) ^cosh{y/qw) J 

and we can take F as 

-1/4 f cosh{y/qw) y/q smh{y/qw) — qw cosh{y/qw)\ 
\smh{y/qw) y/qcosh{^w) — qw sinh{^w) J 

Then /o = EE"^ is a (geodesic) line when |g| = 1, and is a (hyperbolic) cylinder when \q\ ^ 1. 
Also, /i = FF^ gives CMC 1 Enneper cousins in [10] [26] (see Figure 121). These /i and /o 
deform to each other via the first Weingarten family. 

Now, consider half-lines in the domain (the complex plane C) of the Enneper cousins /i 
emanating from the point w = 0. In all but two directions for these rays, the corresponding 
curve on the surface will converge to a single point in the sphere at infinity dM^. This limit 
point can be different for different directions, and in particular the limit will be one certain 
point (resp. one other certain point) for any ray that makes an angle of less than (resp. more 
than) 7r/2 with the positive real axis in the w-plane. However, in the two special directions for 
these rays where w is pure imaginary, the corresponding two curves on the surface converge to 
infinite wrappings of circles in 9IHI'^. These properties are easily checked, because we have an 
explicit form for F, as above. Furthermore, the behavior is the same for any fj:^^ whenever 
t > 0, which is also easily checked. This behavior is similar to the Stokes phenomenon we will 
see in Example 13. 101 

See Figured] for the first Weingarten family of /o, and Figure [3] for a second Weingarten 
family of /q. 

Example 3.9. To make CMC 1 surfaces of revolution /i, the so-called catenoid cousins [TO] 
[31] ■ one can use g — e^^ in ()3.ip . for /i either real or purely imaginary. For the corresponding 
flat surfaces /o, one can obtain the surfaces called hourglasses (resp. snowmen) by choosing 
real (resp. imaginary) values for /i [IT] [21]. Discrete versions of these flat surfaces can be seen 
in Figure [71 and graphics of the smooth surfaces look much the same, but are smooth. 

Example 3.10. We will now consider the holomorphic function g — = z^l'^ ^ which gives 
the hyperbolic Schwarz map /p for the Airy equation -^j^u — wu = as in Equation (|3.6p . 

The value 7 = f corresponds to the choice q = w, and so is of particular interest. We can 
see this correspondence as follows: We take q to be w. This means we have w = [g'Y, and so 
dw = 2g'g"dz. Then, because dw = [1 / {dg / dz))dz , we have {l/{dg/dz))dz = 2g'g"dz, and so 
the original holomorphic function g, as a function of z, would satisfy 

and so 

g = 

and the scalar factor \J can be removed by replacing z with an appropriate constant real 
multiple of z. 

This surface /o has an umbilic point at w = (so the corresponding caustic will blow out 
to infinity at w = 0, see H?], [IH]j US), and it has a "triangle" of singular points with three 
cuspidal edge arcs connected by three swallowtail singularities, and it has 120 degree dihedral 
symmetry. Similar to the case of the CMC 1 Enneper cousins, starting at the center point 
of the surface (w = in the domain of the mapping) and going in any direction but three 
(i.e. following a ray out from w = in the domain), the corresponding curve on the surface 
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will converge to a single point (one of three possible points) in the sphere at infinity 911^ (as 
in Example 13.81 this limit point can be different for different directions). However, in three 
special directions for these rays, the corresponding curves on the surface converge to infinite 
wrappings of circles in OM"^. These special directions are exactly opposite to the directions of 
the swallowtail singularities. This is called Stokes phenomenon, and was explored carefully for 
this example in [27] , [29] . Portions of this /o can be seen in Figure |4| 

The surfaces ' {0 < t < 1) in a second Weingarten family of /o have the same property 
that the singular set is a triangle of three cuspidal edge arcs connected by three swallowtails. 
Then, as t tends to 1, the three swallowtails converge to the origin w — 0, and when t = 1 the 
CMC 1 surface fi has a branch point of order 2 at w = 0. The deformation near the origin can 
be simulated by the map Yg : {u, v) i— > (j/i, ^2; J/s) given by 

yi = —s[u + V ) + 2uv — |w , ?/2 = 2su + u ~ v , 1/3 = —sv + uv , 

as s = 1 — t tends to zero. 

The first and second Weingarten families for /o can be seen in Figures [5] and [Hj 

More generally, using q — in ()3.1I) will give the flat surfaces /o investigated in [2^ . When 

n is neither —1 nor —2, we have g = c„z "+2 for some c„ > 0. When q = w , then g — 21ogz. 

When q = w^^, we have g = ce~^, and c can take on any value. 

Remark 3.11. It is interesting to note that, in Example 13. 101 while the two choices q = w and 
q = w~-^^'^ give different equations ■£;pu = wu and — w^^^^u, the solutions u to the 

second equation are essentially just w-derivatives of the solutions u to the first equation. A 
computation then shows q and q can produce two flat surfaces that are parallel surfaces of each 
other. In terms of the coordinate z, the corresponding statement is that g = z^l'^ and g = z'^l'^ 
can produce two flat surfaces that are parallel to each other, and in fact this also follows from 
Lemma (taking B to be off-diagonal with off-diagonal entries both z). 




Figure 2. The first Weingarten family associated with a hyperbolic cylinder, 
in the Poincare ball model for H^. Here we take g > constant, so that g' = ^Jq 
and g = qw. Half of each of the surfaces is cut away. Note that the surface is a 
CMC 1 Enneper cousin when t = 1. See Example 13. 81 (Here we have included 
grid lines on these smooth surfaces simply to make the graphics more visible.) 
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t = t = l/4 t=l/2 




Figure 3. A second Weingarten family associated with a hyperbohc cyUnder, 
in the Poincare baU model. See Example 13.81 




Figure 4. Images of the smooth flat surface related to the Airy equation 
— wu = 0, in the Poincare ball model. The left figure is the image of the 
unit disc {\w\ < 1}, shown with grid lines for better visibility. The boundary 
of this left figure is the cuspidal edge curve with three swallowtail corners. 
The middle two figures are the image of the ring {1.17 < w < 2.34}, shown 
twice, once with grid lines and once without. The right figure is the image 
of {3.27 < w < 4.09}, this time shown without grid lines because in this case 
visibility is better without them. See Example 13.101 

4. Discrete CMC 1, flat and linear Weingarten surfaces in 

4.1. Known definition for discrete CMC 1 surfaces. We now describe how discrete CMC 
1 surfaces in were defined in ^13j. For this, we first give a light cone model for in 
Minkowski 5-space M'*'^ that is commonly used in Moebius geometry, and that was used in [13] 
in conjunction with quaternions. Let Ti denote the 4-dimensional vector space of quaternions 
with the usual basis 1, i, j and k, and with the usual notion of quaternionic multiplication. 
Points {xi,X2,X3,X4,xo) g M^'^ can then be given by 

(4.1) x=x,(i _°,)+x2(; _y+-3(o \)+-^[\ J)' 

and the R'*'-'^ metric {X,X) is then given by (/ is the identity matrix) 

{X,X) ■ I = -X^ . 

Now let us view the collection of such trace-free 2x2 matrices with imaginary quaternions on 
the diagonal and reals on the off-diagonal as the set of points in R^'^. We can then define 
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Figure 5. The first Weingarten family associated with the Airy equation, in 
the Poincare ball model. See Example 13.101 The bottom two figures are the 
same CMC 1 surface, shown once with grid lines and once without. 



as a 3-dimensional submanifold of the light cone in this way 
(4.2) If = \ x e R"*'! 



V V * 

This H"^ will have constant sectional curvature —1 when given the induced metric from R"*'^. 
Furthermore, for each point X in the light cone, there is at most one value for r e M so that 
rX lies in H^, so we can alternately view as the projectivized light cone. 

In [T3j, in order to define discrete CMC 1 surfaces in 'M? , the discrete version of the Bryant 
equation 

(4.3) F,-Fp^ Fp , det e R , 

V ^ 9q J 9q — 9p 

was used, where g is a discrete holomorphic function with cross ratio factorizing function apq, 
and p and q are adjacent vertices in the domain D C I? oi g. Note that we have assumed 
9q ~ 9p (see Section 12. 3p . The nonzero real parameter A can be chosen freely. The formula 
for the discrete CMC 1 surface / in was then obtained, analogous to the formula (|2.5[) for 
the case of discrete minimal surfaces, by setting (here L'' is the 4-dimensional light cone in 

(") - r . (-f ) e H= c L< c . whe,e (») = (» j) (j 

where r is the appropriate choice of real scalar to place fi^p in H'^ as defined in (|4.2p . In fact, 

r = 2{hai + ibay^ = ~2{ahi + ialy^ e R . 

Because the entries of F are complex, not quaternionic, it follows that ha is purely imaginary 
quaternionic, so the above matrix /i_p does lie in R^^^, and thus in and then also in H"^ with 
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t = (flat) t = 1/10 




t = l/2 t=l (CMC 1) t = l (CMC 1) 




Figure 6. A second Weingarten family associated with the Airy equation, in 
the Poincare ball model. See Example 13.101 The two lower-right figures are 
the same CMC 1 surface, shown once with grid lines and once without. (The 
upper-left surface here is the same as the upper- left surface in Figure [5l) 

appropriate choice of r. Also, we actually have a 1-parameter family of surfaces, due to the 
freedom of choice of A. 

The solution F is defined only up to scalar factors. This is essentially because Equation 
(|4.3p is not symmetric in p and and can be explained as follows: Consider a quadrilateral in 
D with vertices p, q, r and s given counterclockwise about the quadrilateral, and with p as the 
lower left vertex. Then Fq = Fp ■ % with 21 determined by Equation (|4.3p . Then Fr = Fg • 5B, 
again by ((i?^ . Similarly, Fg = Fp ■ € and = • D. Thus we expect 2t<8 = and a 
computation shows that this is indeed the case. However, we also have Fp equal to Fq ■ 21, by 
(|4.3p with the roles of p and q reversed, and it tuns out that 2t2t = (1 — Xapq)I ^ I. Thus, the 
solution F to (|4.3p is only defined projectively, that is, it is only defined up to scalar factors. 

Nevertheless, because we are scaling by a real factor r in Equation (|4.4p anyways, the re- 
sulting discrete CMC 1 surface is still well-defined. Thus we have seen the following lemma. 

Lemma 4.1. Although the solution F of Equation (j4.3p is multi-valued and only defined up to 
scalar factors, the discrete CMC 1 surface fi given in (|4.4p is well-defined. 

In fact, the upcoming Theorem 14.21 also implies that the discrete CMC 1 surface /i is well- 
defined. 

4.2. New formulation for discrete CMC 1 surfaces. Equivalently to the definition given 
in [TBj . there is another way to define fp, which is given in the next theorem. This new form for 
fp is convenient, because it is clearly analogous to the form used in the Bryant representation 
for smooth CMC 1 surfaces in H'^, and it removes the need for the Moebius- geometric R''^^ light 
cone model for H'^. 

Theorem 4.2. The above description (j4.4p for the discrete CMC 1 surface given by F is equal 
to the surface given by p FF'^ , up to a rigid motion ofM^. 
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Proof. The matrices fp in H^, as described in (14. 4p . will be of the form 

f-{AC + BD)j + i{AD - BC) CC + DD 

^\ AA + BB j{AC + BD)-i{AD- BC) 

where r is a nonzero real scalar and 



F ■ 

To have /i_p G H^, we should take 



A B 
C D 

1 



AD-BC 

This means that the coefficient of the i term in the diagonal entries will be simply ±1. So we 
can view the surface as lieing in the 4-dimensional space M!^'^, by simply removing the matrix 
term with scalar xi from Equation ()4.1|) . 

Now, the projection into the Poincare ball model is 



{x2,X3,Xi,Xo) 

(4.5) 



1+Xq 

{Re{-AC - BD),lm{-AC - BD), ^{~AA -BB + CC + DD)) 



AD- BC+ \{AA + BB + CC + DD) 

On the other hand, if we look at 

1 ^ 1 fAA + BB AC + BD\ _ f yo + ys yi + iy2 

AD-BC AD-BC \CA + DB CC + DD J 

and then project to the Poincare ball, we have 

(2/1,2/2,^3) 



2/1 - W2 2/0 - 2/3 



(4.6) 



1 + 2/0 

{Re{AC + BD), lm{AC + BD), \[AA + BB-CC - DD)) 



AD - BC +\{AA + BB + CC + DD) 
The quantities (|4.5p and (|4.6|) are the same, up to a rigid motion of H'^, proving the theorem. □ 
4.3. Discrete fiat surfaces. To make discrete flat surfaces in , we can now take 

(4.7) Ep^Fp- , det £; = det e R , 

like in (|3.2p for the smooth case. We then use the same formula (|4.4p as for discrete CMC 1 
surfaces to define the discrete flat surface, but with Fp replaced by Ep. In light of Theorem 
14. 2[ the following definition is natural: 



Definition 4.3. For E given as in (14. 7p . where F is a solution of Equation (14. 3|) . 
(4.8) /o = -^EE^ 

is a discrete flat surface. 

Furthermore, in light of the behavior of the normal for smooth flat surfaces in Equation 
, it is natural to define the normal at vertices of a discrete flat surface by 

P 

A discrete version of Equation (|3.4p can then be computed, and becomes 



(4.10) E,-Ep 

\9q~9p 

where A is an arbitrary parameter in R \ {0}. 
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Remark 4.4. Like in Lemma l4.ll for F, this E is only defined up to scalar factors. However, 
the discrete flat surface fo is still well-defined. 

Remark 4.5. Changing gp to d ■ gp for d E M+ gives parallel discrete flat surfaces, and this is 
seen as follows: gp d ■ gp implies 

(4.11, E,^Ei:.E,-{yf ^ 

which implies the original surface /o,p — dct^g ^p-^p changes to 
(4.12) f^p = cosh(log d) ■ fo.,p - sinh(log d) ■ Np . 

This surface /q is a parallel surface of fo, and /q is the same as the original surface /q. It also 
follows from ()4.1ip and ()4.12|) that the geodesic at /o,p in the direction Np and the geodesic at 

/ 1 Q \ T 

Jq p in the direction Np = -j^^^Ep K ij^v same. 

We now give one of our main results: 

Theorem 4.6. Discrete flat surfaces in M.^ have concircular quadrilaterals. 

Remark 4.7. Although they have concircular quadrilaterals, discrete flat surfaces are generally 
not discrete isothermic. This is expected, since smooth flat surfaces given as in Section W% are 
generally not isothermically parametrized as well. 

Proof. Let /o — dct"(g) ^^'^ be a discrete flat surface defined on a domain D Cll? , formed from 
a discrete holomorphic function g on D. Let p = (rn,n), q = (m + 1,?^), r = (m + l,n + 1) 
and s = (m, n+1) be the vertices of one quadrilateral in D, and let Ep, Eg, E^ and Eg be the 
respective values of E at those vertices. We choose a cross ratio factorizing function for g, and 
we denote that function's value on the edge from p to q, respectively p to s, as a^g, respectively 
aps. Then 



Eq = EpU , U = \ X 



1 9q- 9p 
1 



9q-9p 



Es = EpV , = I ^ 



1 9s - 9p 
1 



Er = EpUVi - EpVUi , Ui = \ Jp^ 1 



gr-gs 

1 9r- 9q 



1 



Vi = [ Ac, 

gr-gq 

(note that we have UVi = VUi, by the cross ratio identity in Equation (|2.2p ). 

Now the cross ratio for the quadrilateral given by the four surface vertices /o,p, /o,gj fo,r and 
fo,s is 

C = {fo,p — fo,q){fo,q — /o,r) "^(/o,r — /o,s)(/o,s — /o,p) ^ , 

or equivalcntly 

(/o,p - fo,q)ifo,q - fo,r)^^ = ^(/o,;, - fo.s){fo,s - fo.r) ^ • 

(Note that commutativity does not hold for the product of four terms in this C, so it is vital 
that the order of the product be given correctly.) To prove the theorem, it suffices to show that 
C is a real scalar factor times the identity matrix. 

Note that ^ ^ 

and ^ ^ ^ 

= det{Ep)det{uf^^^^ ^ d^/'^'^^^^^^ ' 
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SO we have 

{hp~foM,,-hr)-'=det{U)Ep{I-^-L-UU^^^^^^ 
and likewise 

So one finds that 

dct(t/i)i;p(dct(C/)(C7^)~i - C/)(det(Fi)V'f ^ - )"Vf ^[/"^ = 

= det(T/i)C£:p(det(F)(t/^)-i - F)(det(J7i)J7f ^ - U^y^Uj^^V-'^ . 
Because UVi = VUi, we have 



1 



The determinants of U and V are real, so for example det(C/)(C/"^) ^ = _ _ fp-s? 

\ 9p- 9q 1 

and we get 

' 



det(C/)(C/^ )-'-U= (|.g, - gp\' + Xa^,) 



1 



\9p-9q 

With similar expressions for the other differences we see that 

dct([/i) \9q - gp\^ + Xapq ( f gq-gr 

det(Vi) |5.-5,P + Aap. ){g,-gr 

_ \9s- gp\'^ + Xaps_^^ ( ( .9s -5r 



So the expression for C is 

/ gq-gr- Ss-Sp 



ttps / 1 



M 9M^g_^]^P a VO 1 

for some real factor ji. This concludes the proof. □ 

4.4. Discrete linear Weingarten surfaces of Bryant type. We can now look at discrete 
linear Weingarten surfaces in H'^, similarly to the approach taken in Subsection 13.41 for the 
smooth case. In the discrete case, one takes Ep{t) = Ep ■ Lp, for 



i+tg^gp , I l+fffpffp \ 

i+t"9pSp 1+f'gpg p 

V 1+*9p9p / 

We can then define the discrete linear Weingarten surface of Bryant type to be 

1 



Note that even though ft has been defined using one particular choice for the deformation 
though linear Weingarten surfaces (and as seen before, this choice is not canonical), the re- 
sulting collection of all discrete linear Weingarten surfaces does not depend on the choice of 
deformation. This follows from properties analogous to those for the smooth case in Section 
331 

• If gm.n is a discrete holomorphic function defined on a domain D, then gm,n + a is also 
discrete holomorphic with the same cross ratios, for any choice of complex constant a. 
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• If a cross ratio factorizing function for gm.n is a, then gm.n defined by 



-a 



pq 



(4-13) gq - cjp -- 

9q-9p 

(where pq represents both horizontal and vertical edges) is also a discrete holomorphic 
function that is well defined on D (i.e. gm,n is not multi-valued once an initial condition 
is fixed in the above difference equation (|4.13p ). again with the same cross ratios as 

Qm^n • 

• If Equation (|4.10p holds, then we also have 
and so both E and 

t) 

will produce discrete flat surfaces in H^, and these two surfaces are parallel surfaces 
of each other (see Remark I4.5p . The two resulting linear Weingarten families will be 
different. 

• For any constant matrix 

i'l 

The fimction 

agni,n + b 



-bgm.n + a 

is also discrete holomorphic, with the same cross ratios as g. 
• If Equation (|4.3p holds, then we also have 



pq 



'9q J gq-gp ' 



and so both F and 

a —b 



F ■ 



b a 



will produce the same discrete CMC 1 surface in H'^, but will give different linear 
Weingarten families. 

Remark 4.8. When taking g (resp. g) to be the discrete power function z'^ (resp. —z"') as in 
p.3p and (12. 4p . Equation (|4.13p will hold with Upq = 1 on horizontal edges and apq = — 1 on 
vertical edges if 

7 + 7 = 2 

(see Lemma 18.11 in Appendix [5]) . Thus by the third item above, g and g will produce two 
discrete flat surfaces that are parallel to each other. We saw this same behavior in the smooth 
case as well, see Remark 13.111 

We have the following result, which can be proven similarly to the way Theorem 14.61 was 
proven. 

Theorem 4.9. For any t, the resulting linear Weingarten surface ft has concircular quadri- 
laterals. 

Remark 4.10. Note again that, like in Remark 14. 7[ when t is not 1, the surface will not be 
discrete isothermic in general. 
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4.5. Examples. We now give three discrete examples, in parallel with the previous smooth 
examples 13. 8[ 13.91 and 13.101 The third example is in the next section \5\ 



Example 4.11. Here we discretize Example 13.81 We define, for any nonzero constant g > 0, 

5m,n = q- {m + in) , 

and then 

{gm+l.n — 5m,n)(5m+l,n+l — 3in+l,n) ^{gm,n+l ~ 5m+l,n+l ) (ffni," ~ 5ni,n+l) = ^1 i 

SO we can define the cross ratio factorizing function as (where rip denotes the n coordinate of 
We take a solution F of Equation (14. 3p . We have 

F ^^ ^ F U F ,i=F V 

where 

1 + A(to + in) —Xq{m + in)(l + m + in) 
Xq^^ l-A(l + m + m) 

+ X{im — n) —iXq{m + in){i + m + in)^ 
iXq^^ 1 — X{im — n — 1) 

We have the necessary compatibility condition 

Then, using this i^m.m we can construct the Weingarten family for this holomorphic function 
5m, n- There are special isolated values of the scaling A that give atypical results, just like in 
the smooth case in Example 13. 81 where |(?| = 1 gave the atypical result of a geodesic line for the 
resulting flat "surface". However, usually the resulting flat surface will be a discrete cylinder. 



Um^n 





Figure 7. Discrete flat surfaces of revolution and their caustics, shown in the 
Klein model for tf. The left surface is called a "snowman", and the one on 
the right is called an " hourglass" . 



Example 4.12. Like in the smooth case (Example I3.9p . we can use discrete holomorphic 
exponential functions gm,n = e'^'^™+™^ for any nonzero constant c e R U [iR) to construct 
discrete flat surfaces of revolution. See Figure [7l 

Remark 4.13. In Figure [71 we use the Klein model for H^. We sometimes find the Klein ball 
model to be more convenient than the Poincare ball model, because geodesies in the Klein 
model are the same as Euclidean straight lines. Since we are dealing with discrete surfaces 
with geodesic edges, this can be convenient. This is particularly useful when looking at the 
intersection set of two discrete surfaces. 
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Figure 8. A discrete flat surface in the Klein model constructed from the 
discrete holomorphic function z^/^, and thus related to the Airy equation. The 
surface is on the left, and the surface together with its discrete caustic is on the 
right. Because the discrete version of is constructed from circle packings, 
every fourth vertex has two adjacent horizontal (resp. vertical) edges for which 
Ms — 9q~ 5pl is the same. Since \dg\ determines the distance to focal points, 
we always have three normal geodesies off vertices of the surface meeting at a 
single focal point. It follows that C/ consists of triangles in this case. 

5. An EXAMPLE RELATED TO THE AlRY EQUATION, AND StOKES PHENOMENON 

The following Example l5.1l is of interest because it has similar properties to the corresponding 
surface in the smooth case: the surface in the smooth case has trifold symmetry and has three 
swallowtail singularities connected by three cuspidal edges, and also has a Stokes phenomenon 
in its asymptotic behavior |27] (see Example 13. lOp . 

5.1. The discrete flat surface made vifith the discrete power function z'^/'^. We now 
give an example related to the Airy equation. 

Example 5.1. For this discrete example (see Figure[S]), we need the discrete version of g = z^/^, 
as defined in Section ^17^ Because we have a discrete holomorphic function z^^^, we are in a 
position to be able to consider the resulting discrete flat surface via Equations (14. 3p , (|4.7p and 
(|4.8p . and also a discrete Airy equation, analogous to Equation (|3.5p . 

Numerical evidence (see Figure [5]) suggests similar corresponding swallowtail singularities 
for this discrete flat surface, like for the smooth case in Example 13. 101 

In Figure O portions of the hyperbolic Schwarz image (i.e. this discrete flat surface) related 
to the Airy equation are shown. The image of the discrete half-line z > (z G Z) is drawn in 
red, which has a single limit point in the boundary sphere. The image of the discrete half-line 
iz < Q {z & ill) is drawn in blue, which has no single limit point in the boundary sphere 
and instead wraps around infinitely many times. This behavior is typical for the continuous 
Airy case along the Stokes direction, and the corresponding curves on the smooth flat surface 
associated with the Airy equation behave in the same way (see Example I3.10p . This provides 
a numerical confirmation of a Stokes phenomenon for the discrete Airy function. 

We also note that the numerics suggest a behavior similar to that of the smooth case re- 
garding singularities, as the image looks from a distance as though it has three cuspidal edge 
arcs connecting at three swallowtail singularities. Furthermore the discrete caustic of this flat 
surface has a similar behavior to that of the corresponding smooth caustic, in that it points 
sharply outward at its center of symmetry. (Caustics will be introduced in the next section.) 



6. Caustics of discrete flat surfaces 

In this and the next section (and in Appendix [9]) , since we will consider flat surfaces and 
their caustics exclusively, we will abbreviate the notation "/o" to "/". 
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Figure 9. Asymptotic limit of two curves in the discrete flat surface associated 
with the Airy equation. 



t = (flat) i=f/fO 




t=l/2 t=l (CMCf) 




Figure fO. The linear Weingarten family of discrete surfaces for a discrete 
flat surface associated with the Airy equation, in the Poincare ball model. (The 
upper- left surface here is the same as the left-hand surface in Figure [H) 



6.1. Definition of discrete caustics. Let C/ be the collection of focal points, i.e. the focal 
surface, also called the caustic, of a smooth flat surface / in H'^ (note that we should assume 
/ is a flat front, so that the caustic will exist). If Ef is the lift oi f — EfEj (determined from 
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g), then (see [IT]) 

is a lift of C/ = EcfE]j^. Although we have just described the caustic in terms of Weierstrass 
data, it is independent of the choice of that data. We know that Cf is a flat surface (see [17], 
|18j). because its lift Ec^ satisfies the following equation: 

(6.2) {Ec,r'dEc, ^ ^ jdz. 

For the case of a discrete flat surface / in EII'^ with discrete lift E, we must first consider 
how to define the caustic Cf. We can define the normal Np as in (|4.9p at each vertex fp of /, 
so we have normal geodesies emanating from each vertex, and we can consider when normal 
geodesies of adjacent vertices will intersect. Once we have those intersection points, we will see 
that we can consider them as vertices of C/ , giving us a definition for C / . 

Lemma 6.1. Let f be a discrete flat surface in with lift E, as in (j4.10p . constructed using 
the discrete holomorphic function g. Let apq be a cross ratio factorizing function for g. Then 
the normal geodesies in M.^ emanating from two adjacent vertices fp and fq will intersect if and 
only if Xapq < 0, in which case the intersection point is unique and is equidistant from fp and 
fq. Furthermore, even when the two normal geodesies do not intersect, they still lie in a single 
common geodesic plane. 

Proof. Consider one edge pq. Applying an isometry of H'^ if necessary, we may assume without 
loss of generality that 

Then, with dgpq ^ gq ~ gp, 




Thus 



Eq = Aap, , , det = 1 - AOp 



The condition for the two normal geodesies to intersect is that there exist reals ti and ^2 so 
that 

A 0\ , A \ coshi2 -T sinht2 /l \ 

In other words (we now abbreviate dgpq to dg, and apq to a), 

/coshti + sinhti \ c^fHT,--, i^A 

n 1,^-1.-, ] ^ cosht2 ■ EqE^ + smht2 ■ Eq { ,^ 



coshi,-sinh^J--^^^-^^^^+^i^l^^-^<0 -ij^^ 

= cosht2- i ^ ,2 2 +smhi2- i i# 2 

There exists a ^2 so that this last sum on the right-hand side is a diagonal matrix if and only if 

\a — \dg\'^ 



> 1 



Xa + |dg|2 

So the normal lines intersect in H'^ if and only if 

\Xa-\dg\^\ > \Xa+\dgf' 
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This is equivalent to 

Aa < 

and then ^2 satisfies 



sinh(t2) = ' ' , cosh(t2) - ' ^' 



\/-A\a\dg\ ^/-A\a\dg\ 
Now, to get the diagonal terms in the above matrix equation to match, we want ti such that 



cosh(ti) 
sinh(ti) 



2(1 - Aa) 
1 

2(1 - Aa) 



2 - Msl"" - TTTTTJ ) sinhi2 + ( \dg\^ ~ tj^ ) coshts 



\d9?j ^ V' m 

Such a does exist, and in fact a computation shows that this ti is equal to ^2- 

Now the sign of <2 = ti determines which side of the quadrilateral the intersection lies on, 
and ^2 = ^1 > if and only if \dg\^ + Aa > 0. 

Note that because ti = t2, the distance from the intersection point of the normal lines to 
either of fp and fq is the same. 

By a further isometry of H'^ that preserves fp and Np, we may change g to e^^ g for some 
constant 6* G M. Thus without loss of generality we may assume dgpq € R. It is then clear from 
the above equations that the two geodesies emanating in the normal directions from fp and fq 
both lie in the geodesic plane {x2 = 0} n H'^ of H"^, with represented as in Section [XT] This 
proves the last claim of the lemma. □ 

Remark 6.2. By an argument similar to that of the proof of Lemma |6. 11 but simpler, we also 
have the following statements: Let / be a discrete flat surface in H'^ with lift as in (|4.10p . 
constructed using the discrete holomorphic function g. Let Upq be a cross ratio factorizing 
function for g. Let Np denote the normal as in (|4.9p at fp. Then the lines in R^^'^ (not H'^) 
emanating from two adjacent vertices fp and fq in the directions of Np and Nq (respectively) 
will either be parallel or will intersect at a unique point that is equidistant from fp and fq. The 
distance from either fp or fq to that intersection point is 

\9q - 9p? + ^"p 



\9q-9p? - ^'OLpq 

This is true on each edge pg, regardless of the sign of Xupq. Furthermore, these two lines 
stemming from fp and fq will not be parallel if 

^ j \9q-9v? 

Now, for all that follows, we introduce the following assumption. 

Assumption: For g the discrete holomorphic map used to construct a discrete 
flat surface, assume that all quadrilaterals in the image of g in the complex 
plane are properly embedded. 
This implies the following properties: 

(1) dgpq is never zero for any edge pq between adjacent vertices p and q in the domain D. 
Note that this was already assumed in Section [2731 

(2) dgpq and dgps are never parallel for any square of edge-length 1 with vertices p, q, r, s 
in D. (If dgpq and dgps are parallel, then gp, gq, gr, gs all lie in one line, which implies 
that the interior of the quadrilateral in the complex plane is a half-plane, which is not 
properly embedded.) 

Because the cross ratio for g is always negative, the term Aa^^ is negative on exactly all of 
the horizontal edges in the domain D dl? of /, or exactly all of the vertical edges. Applying 
a 90 degree rotation to the domain if necessary, we can, and now do, assume without loss of 
generality that Xapq is negative if and only if the edge pq is vertical in D. 
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Thus every vertical edge pq provides a unique intersection point of the normal geodesies, 
which we denote by {Cf)pq. The domain Dq of this new mesh is now the collection of all 
vertical edges of D. If we let the vertical edges be represented by their midpoints, then we can 
consider C/ to also be defined on a square grid, but now a shifted one, i.e. Dc CZx (Z+ i). 

Definition 6.3. We call the discrete surface 

Cf — {{Cf)pq \ pq is a vertical edge} 

the caustic, or focal surface, of /. 

This new discrete surface Cf is generally not a discrete flat surface, as the vertices in the 
quadrilaterals of the caustic will generally not be concircular. However, it has a number of 
interesting properties closely related to discrete flat surfaces, which we will describe. The first 
property is stated in this lemma: 

Lemma 6.4. The quadrilaterals of (Cf)pq lie in geodesic planes ofM?. 

Proof. This is clear from the geometry of the construction, and from the fact that adjacent 
normal geodesies always lie in a common geodesic plane of H"^, by Lemma l6.ll □ 

6.2. Discrete extrinsic curvature: first approach. Now that we have seen the proof of 
Lemma 16.11 we are able to give one geometric justification for why we can say the discrete 
surfaces in Section are "flat". We will give an argument here showing the "discrete extrinsic 
curvature" is identically 1. 

For a smooth flat surface given by a frame E solving (|3.4p . we find that the two functions 



give the principle curvatures of the surface, and also give the inverses of the distances to the 
focal points in R^'^ of the curves in the surface along which either y or x is constant. Here we 
are considering the focal points found in M?'^ (not H^), but we are finding those focal points 
with respect to the normal directions to the surface in given by the normal vectors (|3.3[) . 
which are actually tangent vectors to itself. Since the extrinsic curvature fcifc2 is exactly 1, 
if ki ^ki, we have |fci| < 1 and \ki\ > 1, and we have that 

(6.3) arctanh (fci) — arctanh (kn^) = arctanh | — r^rl- I — arctanh | — r^rl- I = . 

This last right-hand equality, of course, is clear from the fact that \gx\ = \gy\- However, the 
right-hand equality encodes that the extrinsic curvature is exactly 1 in a way that can be applied 
to the discrete follows: the corresponding equation in the discrete case is given by the 

corresponding summation about the four edges (assume Xapq — Aa^s > and Xaqr — Xctgp < 
- the other case can be handled similarly) 

(6.4) arctanh ( +, '^^ ' ) - arctanh ( ^^'^'^^ + \gr - 9,\' 



Xapq + \gq - gp\'^ J \-Xaqr + \gr-g 

arctanh _ arctanh ( ^^^p + Ig, ~ 9s\\ , ^ ^ 



Xars + \gs - ffrp J \~Xasp + \gp - gs 

of each quadrilateral with vertices associated to p, q, r, s (given in counterclockwise order about 
the quadrilateral in D) in the discrete surface. The analogous geometric meaning of 

-Xapq + \gq -ffpP 
^"pg + \9q - 9p? 

is preserved in the discrete case, as it is the inverse of the (oriented) distance from either fp or 
fq to the intersection point of the geodesies in R'^'^ stemming off of fp and fq in the directions 
of Np and Nq, respectively (see Remark l6.2p . Furthermore, Equation (|6.4p follows immediately 
from the definition of the cross ratio factorizing function a. In this sense, we can say that the 
"discrete extrinsic curvature" is identically 1. 
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6.3. Discrete extrinsic curvature: second approach. We now consider a second approach 
to discrete extrinsic curvature. For a smooth surface of constant extrinsic curvature 1, the 
infinitesimal ratio of the area of the Gauss map to the area of the surface is exactly 1. So 
another way to give a notion that the "discrete extrinsic curvature" be identically 1 for a 
discrete flat surface is to show the analogous property in the discrete case. That is the purpose 
of the following lemma. 

Lemma 6.5. Let f be a discrete flat surface in with normal map N . Consider the vertices 
fp, fq, fr, fs of one quadrilateral (in H^j of the surface associated with the quadrilateral with 
vertices p, q, r, s (given in counterclockwise order) in D. These four vertices fp, fq, fr, fs also 
determine another quadrilateral T ^ , now in M^'^, again with vertices fp, fq, fr, fs, but now with 
geodesic edges fpfq, fqfr, frfs, fsfp mM'^'^, and which is planar inM.^'^. Likewise, the normals 
Np, Nq, Nr, Ns determine a planar quadrilateral Tn in K^'^ with vertices Np, Nq, Nr, Ng and 
with geodesic edges NpNq, NqNr, TVpVT, NsNp in R^^^ 

These two quadrilaterals Tj and Tn lie in parallel spacelike planes o/M'^'^ and have the same 
area. 

Proof. Because fp, fq, fr, fs he in a circle C in H'^, there exists a planar quadrilateral Tf in M^'^ 
with edges that are geodesies in M?'^ , and with vertices fp, fq, fr, fs- 

Since Np and Nq have reflective symmetry with respect to the edge of J-f from fp and fq (see 
Remark 16. 2p . and since similar symmetry holds on the other three edges of J-f, we know that 
fp + Np, fq + Nq, fr + Nr, fs + Ns are the vertices of a planar quadrilateral J-f+N with geodesic 
edges in M.^'^. It follows that Np, Nq, Nr, Ns are then the vertices of a planar quadrilateral 
with geodesic edges in M.^'^. In fact, J^f, J-f+N and J^n all he in parallel spacelike planes. 

The goal is to show that J^f and J^n have the same area. Since J^f and J^n are parallel, it 
is allowable to replace the metric of R'^'^ with the standard positive-definite Euclidean metric 
for and simply prove that J^f and J^jy have the same area with respect to that metric. The 
advantage of this is that it allows us to use known computational methods involving mixed 
areas. See for example, for an explanation of mixed areas. 

Noting that f ± N lies in the 3-dimensional hght cone of R3'\ i.e. {f + N, f + N) = 
{f ~ N, f — N) — 0, we define the lightlike vectors 

Gi = i(/ + iV), G2 = i(/-iV). 

(In fact, the normal geodesic at each vertex of / in is asymptotic to the lines in the light 
cone determined by Gi and G2.) The two concircular sets {(Gi)p, (Gi)q, (Gi)r, (Gi)s} and 
{(G2)p, {G2)q, (G2)r, {G2)s} havc the same real cross ratio, so the two quadrilaterals in R^'^ 
that they determine are either congruent or dual to each other. (In fact, it is generally the 
second of these two cases that occurs.) It follows that the mixed area of Gi and G2 is zero [5]. 
We then have, with "A" denoting area and "MA" denoting mixed area, 

A{f) = A{Gi + G2) = A{Gi) + A{G2) + 2AfA(Gi, G2) - A{Gi) + AiG^) , 
A{N) ^ A{Gi - G2) = A(Gi) + A{G2) - 2MA{Gi, G2) = A{Gi) + A{G2) . 
Hence A{f) = A{N). □ 

Remark 6.6. The proof of Lemma [Ol shows that for a vertical edge pq (of length 1) of D, the 
geodesic edge fpfq, the geodesic through fp in the direction of Np and the geodesic through fq 
in the direction of Nq form the boundary of a planar equilateral triangle in H^. In particular, 
it follows that fp, fq, f^ and are concircular, for any value of d. One can also show that 
/p: fq: fq and fp are concircular even when is a horizontal edge in D. This shows that the 
quadrilaterals formed by the two points of an edge and the two points of the corresponding edge 
of a parallel flat surface are always concircular. This in turn implies that a quadrilateral of the 
surface and the corresponding quadrilateral on a parallel surface have a total of eight vertices 
all lying on a common sphere - forming a cubical object with concircular sides. This gives (see 
^) a discrete version of a triply orthogonal system, that is, a map from Z"^ or a subdomain of 
li^ to M.^ where all quadrilaterals are concircular. 
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6.4. A formula for the caustic. In Lemma 16.41 we gave one property of caustics that is 
closely related to discrete flat surfaces. Here we give a second such type of property, as seen in 
Theorem 16.71 below. 

The equation for the lift i? of / is 




However, since the formula (|4.8p for the surface has a mitigating scalar factor 1/ detE, we can 
change the equation above so that the potential matrix has determinant one, without changing 
the resulting surface, so let us instead use: 

. _ 1 _ 

^ I- TP '^P ' -^P ~ ^P^P ' 

•^det hp 



1 




pq \ dgj,q 



We also now assume that A is sufficiently close to zero so that 



\\apq\ < 1 



for all edges pq. 

We now define, for each vertical edge pq, 



(6.5) E^cu - («^P + bE,) • I 







where a and b are any choice of nonnegative reals such that a + b — 1 (recall the definition of P 
in (|6.ip ). This is a natural discretization of the Eq^ for the case of smooth surfaces (see (16. ip ). 
where we now must take a weighted average of Ep and Eq, and we allow any choice of weighting 
{a,b). For the smooth case in Equation (|3.4p . the fourth root of the upper right term of E^^dE 
divided by the lower left term gives the -y/g' appearing in Equation (|6.ip . For the discrete case 
in Equation (|4.10p . the fourth root of the upper right term of Ep'^(Eq — Ep) divided by the 
lower left term gives the \fdgpql \J\apq appearing in Equation (|6.5p here. This explains why 
we insert the ^jXa^q factors here. Note that ^Xocpq is not real, because \apq < 0. 

Theorem 6.7. The formula 



for the discrete caustic holds for all vertical edges pq, and this formula does not depend on the 
choice of a and b — 1 ~ a. 



Proof. A computation gives 



\d.gpg\ n \ / \d.gpg\ 



{Ec,)pq{Ec,)pq^ - S ■ Ep \^ y^ r^P =d^-^M "o"" VE^)^^ 

\dgpq\ / \ \dapq\ 



2ab- 



pq 

-T 



The scalar factor S is the only part of {Ecj)pq{Ecj)pq that depends on a and 6, but this 
scalar factor is irrelevant in the formula (|6.6p , so we see independence from the choice of a and 
b. The result now follows from the proof of Lemma [6Tl □ 
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Remark 6.8. However, the choice of normal direction at the vertices of C/ does depend on the 
choice of a and b, as the foUowing equation shows: 

1 ^ ^ T 



(6.7) Npg - {Ecf)pq f Q _-|^ j {Ecf)pq - EpilEp , 




Equation (|6.7p imphes that, for both horizontal and vertical edges, the adjacent normal geodesies 
of the caustic typically do not intersect, for any generic choice of a and b. 

7. Singularities of discrete flat surfaces 

The purpose of the next results is to show that the discrete caustics in Section [6] have 
properties similar to the caustics in the smooth case. 

In what follows, we will regard both the discrete flat surface / and its caustic C/ as discrete 
surfaces that have edges and faces in (not just vertices). Since there is a unique geodesic 
line segment between any two points in H^, the edge between any two adjacent vertices of / 
is uniquely determined. The same is true of C/. Then, since the image of the four vertices of 
any given fundamental quadrilateral in D (resp. in Dc) under / (resp. C/) has image lying 
in a single geodesic plane (see Theorem 14.61 and Lemma 16.41) , the image of the fundamental 
quadrilateral in D (resp. Dc) can be regarded as a quadrilateral in a geodesic plane of H'^ 
bounded by four edges of / (resp. C/). This is the setting for the results given in this section. 

In the case of a smooth flat surface (front) / in H'^ and its smooth flat caustic C/, every 
point in Cf is a point in the singular set of one of the parallel flat surfaces of /. In Lemma 
17.11 we are stating that every point in the discrete caustic C/ of a discrete flat surface / is a 
point in the edge set of some parallel flat surface of /. This, in conjuction with Theorem 17.31 
suggests a natural candidate for the definition of the singular set of a discrete flat surface. The 
proof of Lemma l7. II is immediate from the definitions of parallel surfaces and caustics. 



Lemma 7.1. Let f be a discrete flat surface defined on a domain D C 7? determined by a 
discrete holomorphic function g : D C with properly embedded quadrilaterals, and let C / 
be its caustic. Let P £ 'M? be any point in Cf, so P lies in the quadrilateral T of Cf that 
is determined by two adjacent vertices fp, fq of f and the normal geodesies (which contain 
two opposite edges of T ) in the directions Np, Nq at fp, fq, respectively. (Thus pq will be a 
horizontal edge of D.) P can lie in either the interior of T , or an edge of T , or could be a 
vertex of T . 

Then P lies in the edge fpfq of some parallel surface of f . 

The next proposition will be used in the proof of Theorem 17.31 

Proposition 7.2. Let f be a discrete flat surface with normal N produced from a discrete 
holomorphic function g with properly embedded quadrilaterals. Then for all vertices p, Np is 
not tangent to any quadrilateral of f having vertex fp. 

Proof. Take a quadrilateral with vertices fp, fq, fr and of / so that is a horizontal edge of 
D. We may make all of the assumptions in the proof of Lemma |6.1[ including the assumption 
that dgpq is real. Of course, the normal Np = diag(l,— 1) lies in the hyperplane {x2 = 0} of 
R'^^^ (here we regard points of M'^'^ as Hermitean matrices as in Section |3T|) . and so does the 
point fq, since dgpq G R. Furthermore, pqis a horizontal edge of D, so Xupq > 0, which implies 
that fq does not lie in the geodesic in containing fp and tangent to Np. However, g has 
properly embedded quadrilaterals, so both dgps and + dgps will not lie in M, and thus fs 
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will not lie in the hyperplane {x2 = 0}. It follows that Np will not be parallel to the geodesic 
plane in containing the two geodesies from fp to fq and from fp to fs- □ 

For a smooth flat surface (front) / and its caustic C/, a parallel flat surface to /, including 
/ itself, will meet C/ along its singular set, and that singular set is generally a graph in the 
combinatorial sense (whose edges consist of immersable curves). Furthermore, all vertices of 
that combinatorial graph have valence at least two. For example, cuspidal edges form the 
edges of this graph, and swallowtails give vertices of this graph with valence two. In particular, 
no cuspidal edge can simply stop at some point without continuing on to at least one other 
cuspidal edge (as this would give a vertex of valence one). The following theorem shows that an 
analogous property holds in the discrete case. Note that when we are speaking of the vertices 
of this combinatorial graph in the theorem below, these vertices are not the same as the vertices 
of the discrete flat surface, nor its discrete caustic, in general. 

Theorem 7.3. Let f be a discrete flat surface defined on a domain 13 C determined by a 
discrete holomorphic function g : D <C with properly embedded quadrilaterals, and let Cf be 
its caustic. Let be a parallel surface, and let Sd be the set of all P G as in Lemma \7.1\ 
for that value of d, and for any adjacent endpoints p and q of a horizontal edge of D. Assume 
that no two adjacent vertices of are ever equal, and that the faces of Cf are embedded. 

Then Sd is a graph (in the combinatorial sense) with edges composed of geodesic segments 
lying in the image in of the horizontal edges of D under f^, and with all vertices of Sd 
having valence at least two. 

Remark 7.4. The snowman shown on the right-hand side of Figure[7] (see Example l4.12p provides 
an example to which Theorem 17.31 applies. The Airy example in Section [5] also satisfies the 
conclusion of this theorem (see Figure [5]), although it does not actually satisfy the condition in 
the theorem that the faces of the caustic be embedded. 

Remark 7.5. At least one of the assumptions in Theorem 17.31 that the quadrilaterals of Cf are 
embedded and that no two adjacent vertices of Z'* coincide is necessary. Without them, the 
discrete hourglass, as seen in Figure [7] (see Example I4.12|) . would provide a counterexample to 
the result. However, it is still an open question whether both of those conditions are really 
needed. There are reasons why it is not obvious that we can remove one of those two conditions. 
We explore those reasons in Appendix [^1 

Proof. We must show that all vertices of Sd have valence at least two. There are essentially only 
two situations for which we need to show this, one obvious and one not obvious. The obvious 
case is when an entire edge fpfq lies in Sd, and then the result is clear. The non-obvious case 

that we now describe, where only a part of /p/g lies in Sd, is essentially only one situtation, 
since any other non-obvious situation can be reformulated in terms of the notation given below. 
Without loss of generality, replacing /'^ by / if necessary, we may assume that d = \. 

Let D be a domain in I? containing pi = (0,1), p2 = (0,2), = (0,3), p^ = (0,4), 
qi — (1,1), 52 — (1,2), qz = (1,3), 94 = (1,4), and let 5 be a discrete holomorphic function 
defined on D. Let f = f^ he the resulting discrete flat surface. We have a normal direction 
defined at each vertex of /, which determines normal geodesies £p., £q. at the vertices fp., 
fq^, respectively. Note that ^p. and fg. never intersect (and thus /p. and /g. are never equal), 
although they do lie in the same geodesic plane, and that ip. and Ip^^^ (resp. iq. and ^g.^J 
intersect at a single point that we call Cf^p^ ^j^^ (resp. C/^g^^^J, by Lemma [6.11 The point 
^f,Pi,i+i (resp. Cf^q. .^j) is equidistant from the two vertices fp. and /p^+i (resp. fq. and fqi^^), 
as in Lemma [6. II 

Now the caustic Cf has two quadrilaterals, described here by listing their vertices in order 
about each quadrilateral: 

•^1 — (C/,pi2, C'/,gi2, C/_g23, C/_p23) and — {Cf,p23,Cf^q2^,Cf^q^^,Cf^p^^) . 
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Figure 11. A graphical representation of the argument in the proof of The- 
orem [731 



Let a; be a point in the geodesic edge fp^fq^ so that x also lies in the edge Cf^p^^Cf^q^^. Since no 
two adjacent vertices of / are ever equal, it follows that x lies strictly in the interior of /p2/g2- 
(See the left-hand side of Figure[TT]) Thus, since the face fi of the caustic is embedded, there 
exists a half-open interval X = (y, x] or X = [x, y) contained entirely in the interior of /p2/(?2 so 
that T lies in Ti. 

Thus X can become a vertex of the graph Sd=i- We wish to show that the valence at x is at 
least two. This would mean that the visual representation is more like in the right-hand side 
of Figure [TTl where the quadrilateral of / with vertices /p2 j /(32 > /ga : /ps is nonembedded. It 
suffices to show that there exists a half-open interval I = (y, x] or X = [x, y) contained entirely 
in the interior of the geodesic edge fp^fq^ so that: 

(1) I lies in the face T2 of the caustic, and 

(2) x = x. 

Because x lies in both the geodesic plane determined by fp^ , fq^^ , fq^ , fp^ and the geodesic plane 
determined by fp^ , fq^,C f^qr^,^ , C f^p^^ , and because Proposition 17.21 implies these two geodesic 
planes are not equal, x must also lie in the line determined by fp3,fq3- Then, because x lies 
in the embedded face jFi, and because both Cf^p^^Cf^q^^ and fp^fq^ lie in the geodesic planar 
region between (p^ and iq^ , x must lie in the edge fp^ , fq^ itself. (Note that we have now proven 
that the quadrilateral with vertices fp^ , fq^ , fq^ and fp^ is not embedded, so in fact the visual 
representation must be more like in the right-hand side of Figure [TT]) 

Keeping in mind that x also lies in the edge Cf^p^sjC/.q^^, then since both edges fpsjfq^ 
and Cf^p23,Cf^q23 lie in the plane determined by J-2, and since T2 is embedded, we conclude 
existence of such an interval 2 with the required properties. □ 

The above Lemma [TTT] and Theorem 17.31 suggest that Sd has many of the right properties to 
make it a natural candidate for the singular set of any discrete surface in the parallel family 
of/. 
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Remark 7.6. We have chosen to consider the set Sd in the image f{D) of a discrete flat surface, 
rather than in the domain D itself (as is usually done for the singular set in the smooth 
case), because Sd becomes a collection of connected curves in the image, while in the domain 
D it would jump discontinuously between points in the lower and upper horizontal edges of 
quadrilaterals of D (as we have seen in the above proof). However, one could remedy this by 
inserting vertical lines between those lower and upper edge points in quadrilaterals of D, and 
then consider the set in the domain. 

8. Appendix: the discrete power function 

In Figure I10| we have drawn some of the discrete surfaces in the linear Weingarten family 
associated with the discretization of the Airy equation. For making these graphics, we used the 
discrete holomorphic power function. We explain here how to solve the difference equation for 
determining the discrete power function, as follows: 

_ {gm,n-gm + l,n){gm+l,n + l-gm,n + l) 



(9r,i+l,n-gm+l,n + l)(Sm,,i + l-Sn 



-1 , 



gm+l,n-gm-l,n ffm . ,i + 1 - flm , n - 1 ' 

with the initial conditions 

30,0 = 0, 51,0 = 1, 50,1 = «^ ■ 
Once we know gm.o and 50, n, the full solution is given by solving the first equation for the cross 
ratio. 

We fix n and set 

9m = 9m,n and Gto — gm,n+i ■ 
Then, it is seen from the cross ratio condition that 

^ _ _ (ffm+l — gm){Gm — g-m) — (.9m+l - QmY 

\gm+l — gm) + (.txm — gm) 

If we set 

(^7n — gm and — g^yi^i g^yi , 

then {am} satisfies 

2 

Prrt^m Prn 

O-m+l 



am ~t~ Pm 



We define the recurrence relations: 



^m+l — Pm^m Pm^m: 
Cm+1 — ~^ PmCm 



SO that 

The initial conditions are 



Co = 1 and 60 = 5o,n+i - .go,™ • 
The relation is written in the form 

bm+l \ ^ ( Pm -Pm \ f ^m 
Cm+1 J V 1 Pm J \ Cm 

and the eigenvalues of the 2x2 matrix just above are (1 ± i)pm- 
The difference equation satisfied by Cm is 

Cm+2 - {Pm +Pm+l)Cm+l + '^PmCm = (c > 0) 

where Co = 1 and ci = go,n+i + gi.n — 2.go,ri- Once we have determined {cm}, then 

1 1 ^m 

^m Cm-\-l PmC-m auQ am 

Cm 
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determine 

Then, using Equation (|2.4p . we can determine gm,n for all nonegative m and n. 
We have the following fact, which was also stated in [2J: 

Lemma 8.1. Suppose that gm,n is the discrete holomorphic function solving ()2.3p and ()2.4p 

for one choice of "f, and suppose gm,n is the same, hut with 7 replaced 6?/ 7 = 2 — 7. Then gm,n 
and ~'gm,n satisfy Equation (j4.13p with Upq = 1 (resp. Upq = —1) on horizontal (resp. vertical) 
edges. 

Proof. That (|4.13p holds on the edges (m, 0)(m+l, 0) and (0, n)(0, n+1) can be easily confirmed 
from Equation (|2.4p . Then an induction argument proves the result on all other edges as 
well. □ 

9. Appendix: On a maximum principle for discrete holomorphic functions 

If g{z) is a smooth nonconstant holomorphic function with respect to the usual complex 
coordinate z for C, then log \g\ is a harmonic function, and the maximum principle for harmonic 
functions tells us that log \ g\ cannot have a local finite minimum at an interior point of the 
domain. Thus, if \g\ has a local minimum at an interior point zo, it must be that 5(^0) = 0. 

There are various ways to discretize the notions of holomorphicity and harmonicity. See [3] , 
[S], [H], 0, [2^, [23], [24], [32], to name just a few of the possible references - however, the 
history of this topic goes back much further than just the references mentioned here. These ways 
do provide for discrete versions of the maximum principle. The definition we have chosen here 
for discrete holomorphic functions based on cross ratios, however, does not satisfy a particular 
simple-minded discrete version of the maximum principle, as we can see by the first explicit 
example below. A more sophisticated consideration is needed to produce a proper discrete 
version of the maximum principle, but we do not discuss that here, as the simplest questions 
are what are relevant to Theorem 17.31 

Example 9.1. Set D {{m,n) \ - 1 < m < 1, -1 < n < 2}. Then, with i = set 
50,-1 I - 6i , 50,0 = \ , .90,1 = \ +i 1 90,2 = 3^ + 10« , 
5-1,0 = I - 5\/35 + ^i + 3cos(-|7r) + 3isin(-|7r) , 
51,0 = I + 5V15+ ii + 2cos(-|7r) + 2isin(-|7r) , 
and extend g to all of D so that all cross ratios oi g on D are —1. That is, take g so that all 
Q!(m,Ti)(m+i,n) = ~1 ^ud all = 1. Thcu we have the following two properties: 

(1) Amongst all vertices of D, \g\ has a strict minimum of 1/3 at the interior vertex (0,0). 

(2) Amongst all edges of D between adjacent vertices p and q (both horizontal and vertical), 
\gq — gp\ has a strict minimum of 1 at the interior edge from (0, 0) to (0, 1). 



Figure 12. Caustics C/ which meet for some d at a single point, with g 
taken as in Examples 19. II and 19.21 The caustics are shown in red, and the 
are shown in blue. 
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Taking g as in the previous example, and taking A = 1/100, we can produce the flat surfaces 
and the caustic C/. It turns out that the surface and the caustic C/ intersect at just 
one point if c? w 1/10 is chosen correctly. See the left-hand side of Figure fT2l In this case, two 
adjacent vertices (coming from (0, 0) and (0, 1)) of do coincide. However, the quadrilaterals 
of the caustic C/ are not embedded in this case, so this example does not suffice to show that 
both assumptions at the end of the first paragraph of Theorem 17.31 are truly needed. 

In light of the equations in the proof of Lemma [6.1[ a natural next step toward understanding 
the role of coincidence of vertices of /'^ in Theorem 17.31 (and thus toward understanding which 
assumptions the theorem really needs) is to consider a discrete holomorphic function with the 
properties as in the next explicit example. 

Example 9.2. Taking the same domain D as in the previous example, we now set 
.90,-1 = 1 - |« , 50,0 = , 50,1 = i , 30,2 = 1 + |« , 
3-1.0 = -\/90+ \i + ^cos(- jiyTr) + sin(- jijTr) , 
31,0 = -\/6+ 5? + |cos(-|7r) + |isin(-f7r) , 
and, like in the previous example, we extend g to all of D so that all a.(^m,n)(m+i.n) = ~1 ^-nd 
all a(m,n)(m.n+i) = 1- Then we have the following two properties: 

(1) Amongst all vertical edges of D between adjacent vertices p and g, \gq — gp\ has a strict 
minimum of 1 at the interior edge from (0,0) to (0, 1). 

(2) Amongst any three vertical edges pq = (—1, n)(— 1, n + 1) and pq = (0, n)(0, n + 1) and 
pq — (l,n)(l,n + 1) at the same height in D, \gq — gp\ is (strictly) minimized at the 
central edge (0, n)(0, n + 1). 

Using the function g in Example 19.21 one might hope that the resulting surface would 
show the necessity of the assumption in Theorem 17.31 that the adjacent vertices of do not 
coincide. However, it turns out that the quadrilaterals of C/ are not embedded in this case as 
well. See the right-hand side of Figure [T2l where again A = 1/100 and d (w 1/10) is taken so 
that the two vertices of coming from (0, 0) and (0, 1) coincide. 

Because of these sublties, we leave open the question of whether just one of the two conditions 
in Theorem 17.31 that 

(1) the adjacent vertices of never coincide, and 

(2) the caustic has embedded faces 

would suffice. 
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